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^ Abstract 

I We compare the nature of electromagnetic fields and of gravitational fields in linearized general relativity. We 

carry out this comparison both mathematically and visually, hi particular the "lines of force" visualizations 
W) of electromagnetism are contrasted with the recently introduced tendex/vortex eigenline technique for visualizing 
gravitational fields. Specific solutions, visualizations, and comparisons are given for an oscillating point quadrupole 
^ source. Among the similarities illustrated are the quasistatic nature of the near fields, the transverse l/r nature 
of the far fields, and the interesting intermediate field structures connecting these two limiting forms. Among the 
differences illustrated are the meaning of field line motion, and of the fiow of energy. 
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I. INTRODUCTION 



There are many similarities between electromagnetic (E&M) radiation and gravitational radiation: 
both travel at the speed of light; both carry energy away from their sources; both consist of transverse 
waves with two polarizations. In addition, Einstein's general relativity, the theoretical underpinning 
of gravitational waves, can be put into a form remarkably parallel to Maxwell's electrodynamics, the 
theoretical underpinning of E&M waves. Despite the many similarities, there are important differences, 
and focusing on those differences helps to give a deeper understanding of both kinds of radiation. 

Our mechanism for exploring those differences will be the visualization of the fields. Visualization of 
E&M field lines has proven very helpful to student understanding. The visualization of gravitational 
fields has been a challenge, but the recently developed technique of using "tendex" and "vortex" line^^El 
provides insights that may be of pedagogical value comparable to the plotting of electric or magnetic 
field hues in E&M. 

Both in E&M and in gravitation, the dynamic nature of radiation fields is of central importance to 
visualization, and in both E&M and gravitation, waves carries energy. In E&M it will turn out that a 
definite meaning can be given to the motion of field lines and to the transport of energy. By contrast, in 
the gravitational case, a definite meaning cannot be given to either of these concepts. This will help us 
understand some important ways in which gravitation fundamentally differs from E&M. 



The rest of this paper is organized as follows. In Sec. II A we give a very brief review of E&M theory 
to serve as a basis for comparison with the elements of gravitational theory that we subsequently present 



in Sec. II B Section III develops the principles of visualization of static fields, for E&M in Sec. Ill A and 
for gravitaion in Sec. [Ill B In Sec. Ill, to illustrate both E&M and gravitational static fields, we focus on 
a particular model that will be useful later, in the discussion of radiation. A dipole is the simplest model 
for a source of E&M radiation, but a gravitational dipole cannot generate radiation. We therefore choose 
the simplest configuraton that can generate both E&M and gravitational waves: a point quadupole. 

For simplicity of visualization, we eliminate the issues of visualizing truly three dimensional fields by 
having the point quadrupole be axisymmetric. Time changing sources give rise to radiation fields for both 



E&M and gravitation. These dynamical fields and their visualization are discussed in Sec. IV with the 
example of oscillating E&M and gravitational quadrupoles. We summarize and restate our conclusions 
in Sec. El 

A few words are in order about the choices that have been made for notation and conventions. 
The principle has been to produce a paper that can be understood by a reader with a minimum of 
mathematical preliminaries. To that end we have avoided certain practices that are common in advanced 
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literature. The following points deserve particular notice, (i) Papers involving relativity typically assume 
units in which the c, the speed of light, is taken to be unity. In order to have expressions in which the 
dimensionality of quantities is more transparent, we do not make that choice; all factors of c explicitly 
appear, (ii) It is common to use the "Einstein summation convention," in which summation is assumed 
for any repeated index. Adopting this convention would allow us to drop the explicit summation symbols 
m Eq. ^ and many subsequent equations. We have chosen, however, to have these summation symbols 
appear, (iii) To avoid the mathematical baggage of covariant differentiation we have been explicit in 
using only Cartesian coordinates and Cartesian components where expressions involve differentiation, as 
in Eq. ([5]). (iv) We avoid "coordinate bases" that are commonly used in computations with tensor fields. 



Rather, the components expressed, e.g., in Eq. (20), are with respect to the familiar spherical coordinate 
orthonormal basis, not the coordinate basis. One of the simplifications following from this choice is that 
indices on components are the same whether they are superscripts or subscripts; their location is chosen 
for convenience. 



II. ELECTROMAGNETIC AND GRAVITATIONAL FIELDS: INTRODUCTORY THEORY 
A. Electromagnetic fields 

For comparison with the gravitational case, it is useful for us to mention the very roots of electro- 
magnetic physics. The electric field E and magnetic field B are defined through the expression for the 
acceleration a due to the total electromagnetic force, the Lorentz force, acting on a point particle of mass 
m and charge q, moving at velocity v,EI 

a= — (E + vxB). (1) 

m 

The fields obey the Maxwell equations, which in rationalized MKS units, take the form 

V-E = VxE+ — = (2) 

ot 

1 

VxB-^— = V-B = 0. (3) 

ot 

Here we have simplified the equations by assuming that they apply in a region devoid of sources and 
of material properties, i.e., we take the charge and current density to be zero, and we assume vacuum 
values of the dielectric constant and magnetic permeability. It will be useful to rewrite these equations 
in component form (in which we assume a Cartesian basis): 
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Here the summations are over the indices of the Cartesian components and coordinates, {x^,x'^,x^}, 
or {x,y,z}. The symbol Cijk is the three-dimensional alternating symbol used in the construction of 
determinants and cross products. It vanishes if any of its indices is repeated (e.g., 6221 = 0), equals +1 
for any even permutation of 1,2,3, or x,y,z, (i.e., 6123 = £312 = £231 = 1,), and equals -1 for any odd 
permutation (i.e., £213 = £321 = £132 = !)• 

The simplest solutions of the Maxwell equations are the time independent solutions, especially the 
point multipole (point charge, dipole, quadrupole ...) solutions for electrostatic or magnetostatic fields. 
Of greatest interest in this paper will be the not-so-simple radiation solutions in which time variation 
is essential. Of considerable importance to these radiation solutions is the Poynting vector, the flux of 
electromagnetic power per unit cross sectional area 



In a general treatment, the computation of the radiation produced by a given distribution of time- 
changing charges and currents leads to retarded integrals over those sources. Here, however, we are 
primarily interested in the description and visualization of the resulting fields, so we simply invoke 
radiation fields without being specific about the internal details of their source. For such purposes the 
choice usually made is that of a point dipole, but for our purpose here this choice has the disadvantage that 
it has no gravitational analog. The lowest order multipole for gravitational radiation is the quadrupol^. 
Accordingly, as our example of an electromagnetic radiating source we choose a point quadrupole, and we 
make the description and visualization as simple as possible by taking the quadrupole to be axisymmetric. 

B. Gravitational fields 

To understand what is meant by "gravitation" in relativistic theories it is best to start with the 
simplest case, gravitostatics, gravitation for static configurations. In this case, if the gravitational fields 
are typically weak (if they do not drive particles to speeds comparable to c) then Newtonian ideas 
can be adopted, with minor modification, to relativistic gravitation. By "gravitational field," in this 




E X B . 



(7) 
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approach, we do not mean, e.g., the downward acceleration at 9.8m/sec^ near the surface of the Earth. 
More generally, if is the usual Newtonian potential, V<I>g is not considered "true" gravitation. Since 
it affects all particles identically, its effects on particles undergoing their natural, freely falling motion 
disappear in a freely falling frame, the inertial frame in the relativistic view of spacetime. 

Gravitation, in the relativistic viewpoint, is the way in which the the natural free-fall motions vary 
from place to place and time to time. In a static configuration, one in which there is no change from time 
to time, the information is contained in the way in which V<l>g varies from place to place. The rate of 
variation of a vector is a tensor, in the case of V<I>g it is often called the gravitoelectric fielcP'. (Because 
this tensor describes the raising of tides on astrophysical objects in a Newtonian setting it is also called 
the tidal tensor.) In a Cartesian basis the content of this tensor is the set of tensor component^ 

The trace of this gravitoelectric tensor is a familiar quantity 

k k 

Just as the electric field is divergenceless outside sources, the gravitoelectric field is traceless outside 
sources. 

The definition in Eq. ^ is valid even in the presence of sources. In that case the right hand side 
has the familiar value AttG times mass density. This equation, then, is analogous to Coulomb's law for 
electromagnetism. It gives a definition of the field in terms of its sources, but only for a static field. 
To deal with radiation we need more general definitions, definitions based on the manifestations, the 
physical effects, of these fields. In the case of electromagnetism, this more general definition is given by 
Eq. ([T]) or, equivalently, Eq. Q. In relativistic gravitation there is no concept of force per se. Rather, the 
manifestations of gravity are seen in the effects on two point particles separated by a small displacement 
s, and with a relative velocity v. The relative acceleration of the two particles, given 

^ = - E - 2 E ^jkpBpmv's"' , (10) 

k k,p,m 

defines the tensors £ and B. The similarity to the Lorentz acceleration in Eq. Q is striking, especially if 
one considers SjkS^ and Bj^s^ to be vectors. There are, of course, differences of detail, one of which is very 
fundamental. In Eqs. ([T]), or Q, the factor q/m describes the special features of the particle undergoing 



electromagnetic acceleration. By contrast, in Eq. (10) there is no reference to any characteristic of the 



particles. In accordance with the so-called "equivalence principle," gravitation acts in the same way on 
any particle. 



The gravitoelectric and gravitomagnetic fields defined by Eq. ( |10[ ) are symmetric ( £jk = £kj, 
Bjk = Bkj)- Outside sources tliese fields are traceless (Ylik^kk=Ylk^kk=^ ) and obey the Maxwell- 
like relation^ 

j \ jk jk / 



jk jk 

which can be written as 



dB 

V-5 = Vx5 + — = (13) 

1 f)£ 

V xB- ^— =0 V-B = 0, (14) 

C"' ot 

with appropriate interpretation of the divergence and curl. Note that the divergence can be taken on 
either index (since the "gravito-" tensors are symmetric) and the curl in these equations is symmetrizecP^. 

Two "theoretical" points bear mentioning: (i) Just as the six independent components of E and B 
contain a complete description of the electromagnetic field at a point, the ten independent components 
of the two symmetric traceless tensors £ and B contain a complete description of the gravitational field 
at a point, (ii) The vectors E and B are "gauge invariant." They cannot, for instance, be made to vanish 
by a mathematical choice. In the same sense £ and B are gauge invarianlP^. This is the mathematical 
equivalent of the physical statement that these quantities are directly physically measurable. 



III. VISUALIZATION OF STATIC ELECTROMAGNETIC AND GRAVITATIONAL FIELDS 
A. Electromagnetism 

1. General considerations for electromagnetic visualization 

The electric and magnetic fields are vector fields, and hence in principle are simple to picture. One can 
use small arrows indicating field direction, with arrow length indicating vector magnitude. Alternatively 
one can sketch field lines, curves to which the vectors are tangent, with the density of these lines indicating 
the strength of the field. The limitations of spatial resolution in these methods are well known. Below 
we will be illustrating vector and line fields with a more modern technique: the line integral convolution 
(Lie) method of Cabral and Leedom^^. In this method the brightness or darkness of pixels is correlated 
along field lines. The method produces images with streaks showing the structure of the field lines in an 
intuitively appealing way and with resolution approaching that of the display. 
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2. Static electromagnetic point quadrupole 



As the simplest example of an electromagnetic multipole source that will be generalizable to gravi- 
tation, we choose an axisymmetric quadrupole. A realization of such a source is shown in Fig. [Tj two 
equal positive charges symmetrically arranged on the z axis about a double negative charge at the origin. 
Since there is no net charge, and no favored positive direction, the configuration has neither a monopole 
nor a dipole moment. The static configuration pictured then has a quadrupole as its lowest nonvanishing 
multipole. 




FIG. 1: A simple model of an electric quadrupole. 

In general, the Cartesian components of an electric quadrupole are given bj^^ 

Qij = j (SxiXj - r'^5ij)p{x)(fx , 
where /9(x) is charge density. For our model in Fig. [l| the specific components are 

Qzz — 4c? Q = Qxx — Qyy — i 



(15) 



(16) 



with Qij = for i / j. Note that the quadrupole Qij is itself a tensor describing the charge distribution 
of the source. It is not a tensor field. The electric field produced by that source is the vector field E. 

The model in Fig. [T]is not a "pure" quadrupole; it has multipole moments of order 2, 4, 6. . . . More 
important, it is not a "point source"; it has a characteristic size d. To create a pure point quadrupole we 
use a limiting process analogous to that for defining a point dipole. We shrink d to zero, while keeping 
finite the product Q = d^q. The result is a point source with only a quadrupole moment. 

The static electric field for the point quadrupole is most simply computed from the electrostatic po- 
tential. A formal procedur^^l ^an be used to find the potential directly from the quadrupole components 
in Eq. (16), or a limiting procedure can be applied to the potential of the three point charges in Fig. [Tj 
The resulting electrostatic potential is 

Q 



(3 cos^ e-i) 



(17) 
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FIG. 2: A Lie of the electric field lines of the azimuthally symmetric static electric quadrupole described in the 
text. 

and hence the electric field E = — V$e has spherical components 

= ( - cos2 - M Ee = COS0 sin^ . (18) 

The electric field topology for this case is illustrated in Fig. [T} For this static electric configuration, there 
are no associated magnetic fields. 

B. Gravitation 

1. General considerations for gravitational visualization 

The problem of visualizing the £jk and Bjk fields is a special example of the question: how does one 
visualize tensorial fields. The "gravito-" fields, £jk and Bjk, are tensors, but at least they are the simplest 
nontrivial type of 3-dimensional tensor: they are second rank (two index) symmetric tensors. In this 
sense they are similar to the most familiar tensors of physics: the inertia tensor, the stress tensor, the 
dielectric tensor of an anisotropic material, etc. 

In the case of a tensor like the inertia tensor of an extended massive object, a reasonable visualization 
is the inertia ellipsoid, a three dimensional ellipsoid whose shape shows the directions and size of the 
principal axes of the moment of inertia of the massive object!^. This ellipsoid for a second-rank symmetric 
tensor is very much the analog of the arrow for a vector. 

The moment of inertia ellipsoid describes a single tensor, not a tensor field. While a display of space 
filled with arrows has some usefulness for the visualization of a vector field, the same is probably not 
true of space filled with ellipsoids. But, just as arrows can be connected together to form field lines, 
the principal axes of tensorial ellipsoids can be connected to form a network of lines with visualization 
properties somewhat similar to field lines. The specific technique for visualization of a symmetric tensor 
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Ajk is to find, at each point in space, the principal directions, i.e., the vectors v satisfying the eigenvector 
condition Ajj^v^ = Xv^. If Ajk is symmetric we are guaranteed that three such eigenvectors exist and are 
orthogonal. The three orthogonal sets of lines connecting these eigenvectors, the "eigenlines" for these 
eigenvectors then give us visual information about our tensor field. 

This method has been used for some time in the visualization of stress (a tensor quantity) in fluid 
dynamics and solid mechanicd^^, and has recently been suggested for use in visualizing the gravitoelectric 
and gravitomagnetic field^2Hll_ Those advocating the application to gravitation give the name "vortex" 
lines to the eigenvectors of B due to the role of that field in driving the precession of spin. The eigenvector 
field lines of S are given the name "tendex" lines, suggestive of the role of the tidal distortions associated 
with £. This approach to visualization holds the promise of giving the kind of intuitive insights that may 
suggest what configurations lead to strong emission of power and of linear momentum in gravitational 
waves. 



2. Static gravitational point quadrupole 

We now consider the gravitational equivalent of the configuration in Fig. [T] Here the objects at 
z = ltd are points of mass M, rather than points of charge q. In analogy with Fig. [T| we include a 
negative mass —2M at the centei'^. The quadrupole components for the gravitational case are given 




FIG. 3: A simple model of a gravitational quadrupole. 



by the same integral as that in Eq. (15), but with p now representing mass density. The components 



therefore are those in Eq. (16) with q replaced by M 



The Newtonian gravitational potential for the mass configuration in Fig. [3] is the same as the electro- 
static potential of the charge configuration in Fig. [T] after the replacements g — )■ M and l/Aireo — )• —G. 
The limit of this potential, for d — t- with Q = Md? fixed, is 

GQ 



(3cos2 6l- 1) 



(19) 
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The components of the static gravitoelectric field are straightforward to compute from this potentiaP^. 
For the static gravitational quadrupole, in analogy with the electric case, there are no gravitomagnetic 
field^l^. The components of the gravitoelectric fields, in a spherical basis, are given by 

12GQ , 2 Q ^\ 

trr = 5 — [6 COS 0—1) 

2AGQ . 

CrB = 5 — Sm COS d (zU j 

t00 — c-ee = ^5 sm ti . 

The components S^A and Sgfh are zero by axisymmetry, and the individual components £^ft,,£gg follow 



from the last of Eqs. (20) and from the tracelessness condition Sff,^ + Sgg = —£„■■ 

We now turn to the issue of visualizing these fields. The structure of the components of £ shows 
that two of the eigenvectors will be in the r, 6 plane, and one in the (j) direction. The eigenlines in the (p 
direction are simple azimuthal circles. The much more interesting eigenlines in the rO plane are shown, 
as Lie images, in Fig. [4| 

The comparison of Fig. [4] and Fig. [2] is very instructive and poses a sequence of questions. The 
electrostatic and gravitostatic potentials are identical aside from trivial replacements. Why are the 
visualizations so different? The immediate answer is that in the electromagnetic case we are picturing 
E = — V<I>e) a vectorial quantity with components d^e/dxK In the gravitational case we are picturing not 
the vector field g = — V$g, with components d^/dx^ , but the tensorial quantity £ with the components. 
d'^^g/dx^ dx^ . The component notation correctly suggests that £ is the gradient!^ of g. Why not, then, 
simply visualize gravitational fields with images of g? For that matter why not simply use potentials in 
both cases? 

In electromagnetism we know the answer. The electrostatic potential is useful for electrostatics, but 
the concept does not carry over to time-changing fields, and hence to radiation. For electromagnetic 
radiation the electric field E is a valid and important concept, but it is not simply the gradient of 
a scalar. The same turns out to be true in gravitation. For gravitational radiation £ is a valid and 
important concept, but its components are not the second derivatives of a scalar field. 
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FIG. 4: Line Integral Convolutions showing the eigenvector fields for the vertically oriented static gravitational 
point quadrupole. The top image shows the field for negative eigenvalues, and the bottom image shows the field 
for positive eigenvalues. 

IV. DYNAMIC ELECTROMAGNETIC AND GRAVITATIONAL FIELDS AND THEIR VISU- 
ALIZATION 

A. Electromagnetic fields 

1. General visualization considerations 

While Lie snapshots are very useful, they show the fields only at a moment of time, and do not capture 
the necessarily dynamic nature of the radiation fields. To address this, a new technique. Dynamic Line 
Integral Convolution (DLIC), has been developecP^. Quite aside from the technical challenges to be 
overcome is a fundamental question underlying the notion of field line movement. Animations can show 
how a certain field line changes from one moment to the next, but what do we mean by "a certain field 
line"? How do we relate the field line at one moment to the "same" field line at another moment? How 
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do we put unchanging "tags" on a field lin^211_ n turns out that there is a reasonable criterion for tagging 
field lines. At least in the case of axisymmetric sources, the mathematical nature of Maxwell's equations 
imposes these tags in a natural way. Furthermore, the same tagging can be justified on purely physical 
grounds. 




FIG. 5: Electric flux E through an area whose boundary is moving with velocity v. 

The mathematical property is illustrated in Fig. [5] The flux of electric field through an area A, with 
unit normal n, changes in time according to 



d_ 
dt 



[ E,-ndA= [ ^■ndA+ [ (E x v) • = [ c'^V xB-ndA + [ (E x v) • dl 
J A J A <Jt Jc J A Jc 

= [ (c^B + E X v) -dl. (21) 
Jc 

We adopt flux conservation as our criterion for the motion of electric field lines. This requires that the 
field line velocity v satisfy 

c^B + E X V = . (22) 

For magnetic field lines, the equivalent condition is 

E + V X B = . (23) 

The component of v along the field lines is meaningless, so these conditions give one constraint on the 
components of v perpendicular to the field lines. In the case of axisymmetry, to be considered below, 
that is all we need. 



For a magnetic field configuration, e.g., an oscillating magnetic dipole, we can interpret Eq. (23) as 



the condition that a charged particle experience no net Lorentz force. Particles trapped in tight orbits 
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around magnetic field lines do in fact, move with field lines, so this condition has a very practical meaning 
in many plasma situations: the motion of field lines is equivalent to the motion of electrons trapped on 
tight orbits around field lines. For an electric field these same considerations apply if we replace the 
electric monopole charge motion with the motion of a (hypothetical) magnetic monopole charge. 

In the radiation field, and in many nonradiative configurations, the electric and magnetic fields are 



orthogonal, i.e., E • B = 0. In this case the solution to Eqs. (22) and (23) are, respectively, 

nE X B E X B 

It is interesting that the numerator in both cases is proportional to the Poynting flux, indicating that 
the flux-conserving motion of the field lines is compatible with the transport of electromagnetic energy. 



2. Oscillating electric point quadrupole 



We can convert the static quadrupole of Fig. [T] to a quadrupole source oscillating at frequency u 
by replacing the static distances d by the oscillating distances do + Ad cos ujt. The "amplitude" of the 
quadrupole then changes from Q = qd'^ to 



Q{t) = q{do + Ad cosut)"^ = qdQ + 2qdoAd cos oot + q{Adcoscjjty 



Er 



47reo 



cos {kr — ojt) I 1 



kr 



sin (kr — uot) 



cos 



(25) 



We now take only the part of this expression that oscillates at frequency w, and we define Q = 2qdQAd. 
With this meaning for Q, the spherical components of fields produced by this source ard^S 



(26) 



Qk^ 



47reo r^ 



sin {kr — ujt) I kr 



6 

kr 



+ 3 



COS {kr — ojt) 



cosv sm( 



(27) 



Qk^ 



Aireocr 



1 ) sin {kr — ujt) + cos {kr — cot) 

3^^^z / kr 



cost/ sin( 



(28) 



where k = lo/c. 

In the limit kr ^ the dominant terms in the field are proportional to and these terms agree 



precisely with those in Eqs. (18). This leads to an important insight about visualization: As kr becomes 



much smaller than unity, the dominant terms in the field are the 1/r^ terms, and these terms have the 
form of the static solution modulated by cosut. In other words, at distances from the source much less 
than a wavelength the field is a quasistatic field, a field with the structure of the static field, but with 
an amplitude oscillating in time. 
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Far from the source, in the region kr ^ 1, the story is very different. Here the radial electric field 
falls off as The dominant fields are the 1/r parts of Eg and B^p: 



Ea 



47reo r 



sin {kr — ujt) cos 6 sin 9 



Qk^ 



47reo cr 



sin (kr — cot) cos 6 sin 9 for kr ^ 1 . 



(29) 



These are the radiation fields, the fields in the limit that r is much larger than the radiation wavelength. 
The fact that the radial component does not contribute to this field is simply a statement that the 
radiation fields are transverse, orthogonal to the direction to the source. 





FIG. 6: A Line Integral Convolution snapshot of the electric field given, at time f = by Eqs. (26) 



and (27) for an oscillating electric quadrupolc.The markers show points at one wavelength on either side 



of the origin. At distances from the origin small compared to one wavelength cT the solutions approach 
those of the static solutions in Fig. [2j This is one frame of a complete movie that can be found at 



http://web.mit . edu/viz/gravrad/visualizations/EM/EMintermediate/ 



Figure [6] gives a LIC snapshot of the fields of the oscillating electric quadrupole at a single instant 
of time. The wavelength A = c/cj = cT of the radiation is indicated on the figure, and it is particularly 
interesting to note that the central region of Fig. [6| the region small compared to A, is indistinguishable 
from Fig. [2} In the opposite limit, the regions far from the central region, to the extent that these are 
included in Fig. [6j show the tendency to take the radiation form, with transverse field lines. What is 
most interesting in Fig. |6j however, is neither the near or the far regions, but rather the intermediate 
region. The figure reveals details of the structure of the transition fields that cannot easily be inferred 



from the mathematical expressions in Eqs. (26) and (27). It shows how the field structure makes the 



transition from the very different nature of the near and the far field^. 

A snapshot of a dynamic field, however, necessarily shows only a single phase of the radiation field. 
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A full description requires an animation, and an animation requires a way of identifying field lines at 



different moments of time, as described in Sec. Ill A The identification principle introduced in that 
subsection gives a single constraint on the velocity of field lines orthonognal to the line itself. In the 
axisymmetric case, geometry gives us the other constraint. In particular, the electric field lines in Fig.[l| 



by symmetry, can have no component. The condition in Eq. (22) then completely fixes the line velocity 
V. A DLIC animation of the dynamic electric field lines is available online2A. 



B. Gravitational fields 



1. Oscillating gravitational point quadrupole 



To create a point quadrupole gravitational source oscillating at frequency w, we can perform precisely 
the same procedure as for the electromagnetic case. We take the distances d to have the form do+Ad cos ujt 
and we take the symbol Q to mean 2MdoAd. The nonvanishing spherical components of £ and 13, the 



analogs of Eqs. (26), (27), are found to be 
-4GQP 



£r9 



3 \ , 3 sin (/cr — wt) 

1 + cos {kr-ujt) + - ^ ' 



£d. 



Bo 



-4GQk^ 



-4GQk^ 

cr^ 
-2GQk'^ 



-GQk'^ 



6 



(kry 
3 cos (kr — Lot) 
3 



3 I cos {kr — ojt) 



kr 
6 



kr — — ] sin (kr — cot) 
kr ' 



sm U cos ( 



kr — ] sin (kr — ujt) 
kr 



sm U cos I 



(—3 + /c^r^) cos (kr — cot) 



(kr)'^ I cos {kr — u}t) 



kr 



2kr I sin (kr — iot) 



sm 



kr 



+ 2kr I sin (kr — iot) 



sm 



(30) 
(31) 
(32) 
(33) 
(34) 



As was the case for the dynamical fields in Eqs. (26), (27) there are interesting limits to these expres- 
sions for both kr <^ 1 and kr ^ 1. In the former case, we find that to the leading 1/r^ order £„, £rei 
and £d)S — £gg, have the form of the values in Eqs. (20), modulated by coswt, while the gravitomagnetic 



components vanish to this order. As is the case in electromagnetism, at distances from the source much 
less than a wavelength the fields are those of a quasistatic source, fields with the structure of a static 
field, but oscillating in time. 

In the opposite limit, for A;r S> 1, the dominant components, those that fall off as 1/r, are 

GQu^ 



£<h 



Bo 



cos ikr—ujt) SIT? 6 
r 

GQoj'^ . j.\ ■ 2 a 
cos(fer — WEjsm p. 



(35) 
(36) 
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The form of these radiation fields suggests transverse fields, a suggestion that can be confirmed with the 
mathematics of general relativity and an appropriate set of definitions and constraints. For our purposes 



here it is sufficient to note that the eigenvector directions can immediately be inferred. From Eq. (35), 
and treating £rr as zero, we have that the matrix of gravitoelectric field components is diagonal and 
hence that there are eigenvectors of the gravitoelectric field in the 9 direction, in the (j) direction and 



(for zero eigenvalue) in the r direction. From Eq. (36) we infer that in the radiation zone there are 
two eigenvectors for the gravitomagnetic field in the 6(j) plane, each at 45 degrees from the and the cj) 
directions. Again, there is an eigenvector field, for zero eigenvalue, in the radial direction. 

Snapshots of the eigenlines of £ are presented in Fig.[7j As in the electromagnetic case, these snapshots 
show what we already learned from the mathematics. Close to the source the eigenline field is quasistatic, 
and far from the source it is transverse. Also as in the electromagnetic case, these figures can clarify 
what the mathematics cannot: the form of the fields in the intermediate zone, at distances from the 
source comparable to the wavelength A = cT. In this zone the fields must make a transition from the 
quasistatic small-r form to the radiating large-r form. 

Both in the mathematics and in Fig. [7] we see that the form of the fields in the intermediate zone is 
even richer for the oscillating gravitational quadrupole than for the oscillating electric quadrupole. In the 
latter case the radiation fields fall off as 1/r, while the near-source quasistatic field has the character 
For the gravitational quadrupole the radiation fields have the same 1/r character, but the near-source 
quasistatic fields behave as 1/r^. 

There is an interesting feature of the gravitational fields that is associated with the mathematical 



expressions more than the visualizations. The expressions in Eqs. (26), (27), and (28), for Er, Eq, B^p, 
are identical to those for £„,£re^^r<i) when the change l/47reor^ — )• —2GQ/r^ is made; this is not true 
only in the near-source zone or the radiation zone, but for all values of kr\ This means that, aside 
from a change in constants and a single factor of r, the electromagnetic solution, including radiation 
fields, is completely contained in the nonradiative part of the gravitational solution. Stated in the 
other direction: the electromagnetic solution contains all of the gravitational solution except for those 



components £^(j, — £eg, and Bgrp, that carry radiation (cf. Eq. (35) and (36)). 

This, of course, is not a coincidence, nor is it an idiosyncrasy of the axially symmetric point quadrupole. 
Rather it is a consequence of the structure of classical theories that describe massless fields, i.e., fields 
that propagate at the speed of light. This structure is most apparent when the fields are described with 
the appropriate mathematics, a set of scalars that result from projecting the fields onto a set of basis 
vectors best suited to the analysis of propagating fieldPl. 

This relationship of the electromagnetic and linearized gravity fields gives us an additional possibility 
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FIG. 7: Lie snapshots of the eigenhnes of the gravitoelectric field S for an osciUating point gravitational 



quadrupole described by Eqs. (30), (31 1 and (33) at time t = 0. The markers show the points at a dis- 



tance of one wavelength on either side of the origin. The top LIC shows the eigenlines for one family23 
of eigenvectors; the LIC on the bottom shows the other family. At distances from the origin small com- 
pared to the wavelength cT, the fields, and hence the eigenlines, approach those of the two families of 
eigenlines for the static solution shown in Fig. |4] The top (bottom) image is one frame of a complete 
movie that can be found at |http : //web .mit . ed.u/viz /gravr ad/visualizations/GRAVinter/EposgravNoBkInt/| 



(http: //web.mit . edu/viz/gravrad/visualizations/GRAVinter/EneggravNoBkInt/ ). 



for visualization. If we project the gravitational tensors with the unit r vector (equivalently take the dot 
product of S and IS with the unit radial vector f ) we get two vector fields, one with r, 0, ^ components 
Ern^rdi ^t^i the other with components Brr,Bre, Brcf,- The E field differs from the vector with components 
Err , £re , £r4> Only by a factor of r, and hence the direction of the two fields is the same at any point. The 
same is true for the corresponding B field. The LICs of the E and B, as in e.g., Fig. [6j fields can 
therefore be considered to be LICs of all parts of the gravitational field except the parts transverse to 
the r direction. It is those transverse parts, of course, that carry the radiation. 
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2. General visualization considerations 



In linearized gravity, there is no equivalent of the principles that define or constrain field line motion 
in electromagnetism. What this means is that we can present a sequence of eigenlines, but we cannot say 
which line at one moment corresponds to which line at another. This fundamental problem turns out not 
to be a barrier if we only want to show the qualitative nature of the field pattern motions in gravitation. 
Intuitively useful flow fields can be guessed at qualitatively due to the existence of singularities. 

In plots of lines tangent to vector fields, the singularities are points at which the vector has zero 
magnitude, so that the direction of a tangent line is undefined. In Fig. [6j such points can be seen both 
in the equatorial plane and along the symmetry axis. If the vector field illustrated were the velocity of 
a fiuid, these singular points would be called stagnation points. 

It is simple to argue that lines of the electric field E in electromagnetic quadrupole radiation moving 



according to the rules of Eq. (22) do not cross singularities. Similary, in the gravitational case, singu- 
larities are easily identified, and the motion of the singularities - their displacement from one snapshot 
to another - therefore gives us a coarse visual sketch of the motion of the field lines; once the position 
of the singularities is known at a new time, the remainder of the field lines can be approximately drawn 
in. Different choices of the details of how they are drawn make no difference in the qualitative content 
of the images. 

Note that the eigenline fields of gravitational radiation shown in Fig. [7] are richer in singularities 
than the electromagnetic vector fields shown in Fig. [6| In addition to "stagnation points," the eigenlines 
contain line singularities in the equatorial plane. In sequences of LIC snapshots the motion of these 
singularities, and the understanding that they "drag" the field lines, give approximate meaning to the 
evolution of the eigenlines. Such qualitative visual fiow fields could be used in making the movies, i.e., 
sequences of LIC snapshots, that have been made available onlin4^^. However, we chose not to do that 
even qualitatively, and have taken the underlying fiow field to be zero in those movies. The eye of the 
observer provides a concept of motion in any case, and we felt that it would be overinterpretation to try 
to augment that perception by imposing a qualitative fiow pattern that matches it. Therefore we do not 
evolve the underlying pattern at all in the gravitational radiations movies online. 

What may be the most interesting point about these visualizations of eigenline dynamics is that, 
unlike motion of fiuid velocity fields and electric field lines, there is no "correct" and "true" motion 
of the eigenlines. There is neither a physical process (e.g., the spiraling of charged particles around 
magnetic field lines, the identification of specific fiuid elements in velocity flows), nor a mathematical 
criterion (e.g., flux conservation) available as a basis for deflning or constraining the motion of lines. 
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Exploiting singularity location is the best we can do. 

The intuitive feeling that there should be a well defined meaning to the motion of field lines may be 
rooted in the relationship of electromagnetic field line motion and energy flow, a relationship, expressed 



in Eq. (24), limited to regions, such as the radiation region, in which the E and B fields are orthogonal. 
This suggests the question of whether energy flow could be used as a guide to the motion of eigenlines, 
and hence the question whether there is a gravitational analog of the Poynting flux. 

There is, in fact, a pragmatic quantitative measure of energy flow in gravitational waves, the Landau- 
Lifschitz pseudotensoi^, but this measure is neither definitive nor useful for our purposes. The fact that 
it is not definitive is important. Just as gravitational acceleration near the Earth surface vanishes in 
a freely falling elevator, many other aspects of gravitation vanish in an appropriately chosen reference 
frame. A consequence of this is the fundamental impossibility of localizing energy in a gravitational wave. 
The Landau-Lifschitz pseudotensor can give only a sort of average energy flux over several wavelengths. 

This already suggests that such a measure is not useful for us; it cannot help us give meaning to line 
motion. The mathematical details confirm that suggestion; the pseudotensor cannot be inferred from 
the gravitoelectric and gravitomagntic fields. In a very rough sense the pseuodensor is constructed from 
mathematical objects that are spacetime integrals of S and T^^. 



V. CONCLUSIONS 



We have compared the mathematics and visualizations of electromagnetic and gravitational fields, 
working up from static field configurations to radiating oscillatory versions of these configurations. Most 
of the details of the mathematics and the images were illustrated with the examples of electric and 
gravitational point quadrupoles. Visualizations, using LICs, were based on the familiar field lines for the 
electromagnetic case. Visualizations of gravitational fields used the recently introduced tendex/vortex 
eigenline formalisnP^. 

We have found, as foretold in the Introduction, that this comparison shows both instructive similarities 
and instructive differences. An important similarity of radiation examples was the change, in both cases, 
from a quasistatic field structure at distances from the point source small compared to a wavelength, 
to a transverse, 1/r, radiation field at distances large compared to a wavelength. The visualizations, in 
both cases, show the field structures that are not easily seen in the mathematics, and show how the fields 
make the transition from the near-source structure to the very different radiation structure. 

An important difference between the two cases lay in the visualization of dynamical fields. The 
motion of electromagnetic field lines has both mathematical and physical meaning, while the motion of 
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gravitational eigenlines, has neither. This difference can be partially ascribed to the difference between 
the nature of eigenlines for a tensor field and "lines of force" for a vector field. But the difference, 
especially regarding energy flow, underscores fundamental differences between electromagnetism and 
relativistic gravitation. 



In Sec. IV B there was a rediscovery and illustration of an interesting mathematical relationship 
between the electromagnetic and gravitational fields. Aside from trivial replacements, the complete 
details of electromagnetic solution - including the radiation fields - is contained within the gravitational 
field solution. To go from the electromagnetic solution to the gravitational, requires only adding the 
gravitational radiation fields. 
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